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Abstract 

The LUCIAE model predictions are compared with experimental data on light vector 
meson {p,uj,(j)) production in pp, p+Pb and Pb+Pb collisions at 158 A GeV. The 
model reproduces the general trends of the data, but fails in describing the yield of 
high-pj^ vector mesons in the central Pb+Pb collisions. Predictions for In+In and p+A 
(A varying from Be to U) collisions at 158 A GeV are also presented. 



1 Introduction 



The enhancement of the strange quark production in heavy ion colhsions, in particular, 
the enhanced yield of meson relative to non-strange p and uj mesons is suggested as a 
possible signature of quark-gluon plasma (QGP) formation [HE]- On the other hand, such 
an enhancement is also predicted by conventional (no QGP) models, in particular, by the 
LUCIAE 0, an updated version of the FRITIOF model HE]. 

The first and most detailed experimental data on the light vector meson (p, tu, 0) production 
in nuclear collisions are obtained at the CERN SPS. However, up to now no attempt was 
undertaken to describe simultaneously these data in the framework of a unique approach. 
In this Note, the LUCIAE model is employed (with somewhat updated parameters) to 
compare the theoretical predictions with the existing data on cf) and (p + uj) production in 
pp, p+Pb and Pb+Pb collisions at 158 A GeV [HI El- In Section 2, the parameters governing 
the strange particle yield and the transverse momentum distribution of produced hadrons are 
described. Section 3 is devoted to the comparison of the model predictions with experimental 
data. Predictions for In+In and p+A (A varying from Be to U) collisions are also presented. 
The results are summarized in Section 4. 



2 The model parameters 

The choice of the LUCIAE model parameters governing the multiplicity and transverse mo- 
mentum distributions of final particles produced in nucleon-nucleon (hadron-nucleon) colli- 
sions at SPS energies is motivated in |H] from the comparison with experimental data. The 
following values of parameters, some of which differ from the default ones [5], will be used 
below: 

- The minimum transverse momentum qTmin transferred during the (semi) hard scattering of 
partons of the colhding hadrons is qmin = VFR(12) = 0.6 GeV/c; 

- The width of the Gaussian distribution of the soft transverse momentum transfer Qt be- 
tween the two colliding hadrons (strings) is given by < Q"^ >= VFR(6) = 0.1 (GeV/c)^; 

- The width of the Gaussian distribution of the primordial transverse momentum Q2T carried 
by the string ends (valence quark or diquark) is given by < Q^t >= VFR(7) = 0.1 (GeV/c)^; 

- The width of the Gaussian distribution of the transverse momentum (p^. andp^ in the string 
reference frame) acquired by hadrons as a results of the string fragmentation is = ay = 
PARJ(21) = 0.405 GeV/c, equal to the default value 0. 

In the LUCIAE model, the broadening of the transverse momentum pt distribution in nuclear 
collisions is suggested to be caused by a novel ('firecracker') mechanism for gluon emission 
by a high-density collective state formed by several overlapping strings j9|. Besides, a phe- 
nomenological mechanism is introduced for the enhancement of the strange quark yield in 
the fragmentation of the collective string states [TU] . 

The model basic parameters governing the yield of hadrons with strange content are: PARJ(2) 
= P(s)/P(u), the rate of the (ss) production compared to (uu) or (dd) production, and 
PARJ(3) = (P(us)/P(ud))/(P(s)/P(d)), the extra suppression of the strange diquark pro- 
duction compared to the normal suppression of strange quarks. 

For 200 GeV pp collisions, these parameters are tuned in by comparison with data on the 
strange particle production: PARJ(2) = 0.2, PARJ(3) = 0.3. At these values, however, the 
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model prediction JT] for the (f) meson average multiplicity < n,^ > in 158 GeV pp collisions 
somewhat (by 18%) overestimates the experimental value ||6j. Below (in Section 3) we will, 
therefore, use for 158 GeV pp collisions a slightly smaller value for PARJ(2) = 0.19, keeping 
the value of PARJ(3) = 0.3. 

Below two options of the model are considered: without and with the final state interac- 
tions(FSI).The latter includes the following inelastic reactions involving p and/or lj mesons: 
nN^pN, NN puj and YN K*uj (F = A or S), leading to an increased yield of cj, 
especially in heavy ion collisions (however, practically not influencing the yield of p), while 
the FSI effects on the $ meson production are neglected [3[TT]. To trace the predicted FSI 
effects concerning the yield of p and uj , the predictions are presented below for both model 
options. 

Note finally, that in the model predictions the following branching ratios for 0, p and uo 
dimuon decay are used: = (2.87 ± 0.22) • 10-^ E^^^ = (4.60 ± 0.28) • 10"^ and 

J^^Hl ~ j^ee _ (j Q'j -|- 0.19) . io~^ The quoted errors of the branching ratios are taken 
into account in calculations. 

3 The model predictions and comparison with experi- 
mental data 

In Table 1 the model predictions for the total average multiplicity < > and the 
< > / < n.,^ > ratio are compared with the NA49 data [H]. An agreement is seen both for 
pp, p+Pb and central (at h < 3.5 fm) Pb+Pb data. The predictions for 158 A GeV In+In 
collisions are also presented. 



Reaction 




<n^> / <n^> 


PP 

NA49 


(1.20 + 0.15) • 10-2 


(0.42 + 0.05) • 10-2 


LUCIAE 


(1.39 + 0.04) • 10-2 


(0.46 + 0.02) • 10-2 


p+Pb 






NA49 




(0.70 + 0.09) • 10-2 


LUCIAE 


(4.7 + 0.2) ■ 10-2 


(0.75 + 0.03) • 10-2 


In+In 






LUCIAE 


0.85 + 0.01 


(0.94 + 0.01) • 10-2 


Pb+Pb 






(central) 
NA49 


7.6 + 1.1 


(1.24 + 0.18) • 10-2 


LUCIAE 


7.00 + 0.04 


(1.13 + 0.01) ■ 10-2 



Table 1: The multiplicity < > and the ratio < n^> j < n^, >= 2 < n^j, > j < n^++?T,^- > 

The model approximately reproduces the pp data 6J on the meson rapidity (Fig. 1) 
and transverse mass (Fig. 2) distributions. The agreement with data for the both 
distributions is somewhat better than in the original paper jTTj, where the predictions for 
the 0-production are presented for the first time in the framework of the LUCIAE model, 
with somewhat different set of parameters mentioned above in Section 2. 
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The comparison with data for Pb+Pb central colhsions [HI is also shown in Figs. 1 and 2. 
The model overestimates by ~20% the data at mid-rapidity and predicts somewhat narrower 
y-distribution (with width ay = 0.97 ± 0.02), than the experimental one (cr^^^ = 1.22 ± 
0.16). The inverse slope of the predicted Mt- distribution, parameterized as dn^/dM^dy ~ 
exp{—MT/Tif)), = 232± 2 MeV, is significantly smaller than the measured one T^^^ = 
305± 15 MeV. The recent attempt [H] to reproduce the Mr- distribution is also unsuccessful: 
the LUCIAE predictions prevail the data at small Mt — < 0.15 GeV/c^, while at large 
Mt — > 0.5 GeV/c^ the predicted values are lower than the experimental ones. 
This inconsistency at large Mt > 1.5 GeV/c^ is more significant when comparing the model 
predictions with the NA50 data for the dimuon channel of production in Pb+Pb collisions 
[?]. As seen from Fig. 3 (the left panel), the model predictions are 1.5-2 times lower of the 
data on the differential cross section integrated over the almost whole domain of impact 
parameter (except for the most peripheral collisions). The predicted value of = 220± 7 
MeV is, however, consistent with the measured one T^^^ = 228± 10 MeV. Better agreement 
is obtained for the p + u: differential cross section (Fig. 3, the right panel). Note, that the 
predicted values of the cross-sections (and of the multiplicities < and < rip+i^ >^^, 

see below) are multiplied by a factor Aqs = 0.5, equal to the NA50 setup acceptance for 
the Collins-Soper angle. In this acceptance, the predicted inclusive cross section (Jp+uj {^^ 
1.5 < Mt < 3.5 GeV/c^ and < i/cm < 1) is equal to 191 ± 6/i6 for the model option 
without the FSI, while, owing to the introduction of the latter, this value reaches 240 ± 8/x6, 
consistent with the measured one. 

The dependence of the vector meson yield on the number Npart of participant nucleons in 
Pb-|-Pb collisions and that per participant nucleon, < /Npart and < rap+^ /Npart, 

are plotted in Fig. 4. It is seen, that the model predicts an increasing yield per participant 
nucleon with increasing Npart, owing to the collective string mechanism incorporated into the 
model. However, this enhancement for (p is much moderate than observed in the experiment 
[7]: the model predictions are lower of the data by a factor about 3 for the most central 
collisions (Fig. 4, the top panels). It should be stressed, that the model fails not only 
for the hidden strangeness sector, but also for the dimuon channel of the high-M^ p + 
uj production. As seen from Fig. 4 (the bottom panels), the model predicts about 1.8 
times smaller < Up^^^ ^fifi^ than measured in the most central collisions, while due to the 
introduction of the FSI this discrepancy reduces to about 1.4. 

The model predictions concerning the 0, p and u production in 158 A GeV In+In collisions 
for the acceptance of the NA50/60 setup (0 < ycm < 1 and Acs = 0.5) are presented in 
Figs. 5 and 6. In this acceptance, the integrated over Mt cross sections are predicted to be 
cT^''=0.21mb, a^''=0.57mb and a^''= 0.83mb for the option without FSI, while (T^''=0.55mb 
and cT^^= 1.09mb when the FSI is taken into account. It is seen from Fig. 5, that the Npart 
- dependence of < > and < n^^ > (when the FSI is taken into account for the latter) is 
significantly stronger, than that for < rii^ > (when the FSI is not taken into account) and 
for < Up >. The inverse slope parameter T (plotted in Fig. 6) for is practically the same 
as for p and uj, when the FSI is taken into account for the latters. 

The predictions for p+A collisions (A = Be, Al, Cu, In, W, Pb, U) are presented, at 
< Vcm < 1 and Acs = 0.5, in Figs. 7 and 8. Again, owing to the collective string 
mechanism, the predicted mean multiplicity < > increases with A faster, than < Up > 
and < >. The ratio < > / < > for the heaviest target (U) exceeds by about 1.5 
times that for pp collisions. The parameter A tends to increase with A (Fig. 8), being for cf) 
systematically higher as compared to that for p and uo. 

?> 



4 Summary 



Predictions of the LUCIAE model, an updated version of the FRITIOF event generator, 
for the p, a;, production in proton-proton, proton-nucleus and nucleus-nucleus collisions at 
158 A GeV are presented and compared with the available experimental data. The model 
parameters used for pp collisions were preliminary tuned by comparison with the data on the 
multiplicity distribution and inclusive spectra of charged particles and the strange particle 
yield in proton-proton (hadron-proton) interactions at the SPS energies. 
The model predictions agree with the data P on the total average multiplicity < > 
and the ratio of < n<^ > / < > for pp, p-l-Pb and central Pb-|-Pb collisions, as well as 
approximately reproduces the Mt and y - distributions of meson in pp collisions. However, 
for central Pb+Pb collisions, the model predicts more narrow y - spectrum and less steep 
Mt - distribution: the model overestimates by 30-50% the yield at low My — < 0.15 
GeV/c^ and underestimates by 20-40% that at higher 0.6 < Mt - < 1.1 GeV/c^. The 
predicted value of the inverse slope parameter = 232 ± 2 MeV of the Mt - distribution 
turns out to be much smaller than the experimental one, = 305 ± 15 MeV. 
On the other hand, the model fits well the shape of the Mt - distribution for the dimuon 
channel of the and p + u production in non-peripheral Pb+Pb collisions [7] . In the high 
Mt - region (1.5 < Mt < 3.2 GeV/c^) the predicted values of = 220 ± 7 MeV and 
Tp+uj = 220 ± 3 MeV agree with experimental ones. The (p + u) yield is also reproduced by 
the model (when the FSI is included into the model), unlike that of meson for which the 
predictions are 1.5-2 times lower than the data. 

The situation worsens dramatically for central Pb+Pb collisions (for the high Mt region). 
The model badly underestimates not only the yield (by a factor ~3), but also that for 
(p + a;) (by a factor ~2). This discrepancy is less significant when comparing with the data 
[0] which concern the — > K^K~ channel in the same range of Mt- 

One can conclude, that the model fits, irrespective to Npart, the total yield of (for pp, 
p+Pb, Pb+Pb), but fails at high Mt region {Mt > 1.5 GeV/c^), especially for the central 
Pb+Pb collisions for which, however, the magnitude of the discrepancy with data is rather 
different for different experiments 

In conclusion, the LUCIAE model reproduces some features of the light vector meson pro- 
duction in proton-proton, proton-nucleus and nucleus-nucleus collisions at the SPS energies, 
but fails in describing the whole totality of the available data. Although the mechanism 
of the formation of multistring states in nuclear collisions, incorporated into the LUCIAE 
model, provides a strangeness enhancement, it turns out to be far from being sufficient to fit 
the data on the high- My meson yield in the central heavy ion collisions. The final state 
interactions, introduced into the model, result in a significant increase of the lj yield and 
lead to a better agreement with the experimental data. 
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Figure 1: Rapidity distributions of (f) mesons in pp and Pb+Pb collisions. The curves are 
the result of the Gaussian fit of the simulated data. 
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Figure 2: Transverse mass distribution of (f) mesons in pp and central Pb+Pb collisions. 
The curves are the result of the exponential fit of the simulated data. 
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Figure 3: Transverse mass distributions of (j) and p + a; in Pb+Pb collisions. The curves are 
the result of the fit(see text) of the simulated data without (sohd lines) and with(dashed 
line) the FSI. 
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Figure 4: Multiplicities < and < Up+^j (left panels) and multiplicities per 

participant nucleon (right panels) in Pb+Pb collisions. 
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Figure 5: Transverse mass distributions of (J), p and lv in In+In collisions. The curves are the 
result of the fit (see text) of the simulated data without (sohd lines) and with(dashed hne) 
the FSI. 
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Figure 6: Mean multiplicities < >^fji < and < integrated over the whole 

Mt - region (left panels) and those per participant nucleon (right panels) in In+In collisions. 
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Figure 7: The A-dependence of the mean muhiphcities < np >, < n^^ > and < 
integrated over the whole My - region and the ratio < n^^ > / < > in p+A colhsions. 
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Figure 8: The A-depcndence of the inverse slope parameters < Tp >, < > and < > 
of the Mt - distributions in p+A colhsions. 
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